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Abstract 

The addition of an adjoint Polyakov loop term to the action of a pure gauge theory at finite 
temperature leads to new phases of SU{N) gauge theories. For SU{3), a new phase is found which 
breaks Z{3) symmetry in a novel way; for SU{4:), the new phase exhibits spontaneous symmetry 
breaking of Z{4:) to Z{2), representing a partially confined phase in which quarks are confined, but 
diquarks are not. The overall phase structure and thermodynamics is consistent with a theoretical 
model of the effective potential for the Polyakov loop based on perturbation theory. 

PACS numbers: 12.38.Gc,11.15.Ha,11.10.Wx 
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INTRODUCTION 



It is well established that SU{N) gauge theories in 3 + 1 dimensions have a low- 
temperature phase in which quarks are confined, and a high temperature phase where quarks 
are deconfined, often referred to as the quark-gluon plasma phase. The deconfinement phase 
transition in pure gauge theories, i.e., without quarks, is understood theoretically as a tran- 
sition between a low-temerature phase where a global Z{N) symmetry is unbroken to a 
high-temperature phase where Z{N) symmetry is spontaneously broken Simulations 
indicate that the transition from confined phase to deconfined phase is similar for all > 3. 
The global Z{N) symmetry appears to always break completely, with no residual unbroken 
subgroup. 

The addition of a term of the form 

- J £x Ha TrAP{x) = -T dt J d^x Ha TrAP{x) (1) 

to the Euclidean action of pure SU{N) gauge theories at finite temperature leads to new 
phases with novel properties. Here P{x) is the Polyakov loop at the spatial point x, given by 
the usual path-ordered exponential of the temporal component of the gauge field Aq in the 
Euclidean time direction. The temporal origin of P is irrelevant due to the trace; because 
the trace is in the adjoint representaton, this additional term respects Z{N) symmetry. 
Of course, this additional term in the action is neither local nor renormalizable in 3 + 1 
dimensions. Thus we must regard this model as an effective theory defined at fixed lattice 
spacing or by some other cut-off. There will be a finite renormalization of the parameter Ha 
in comparing lattice results with continuum. 

This additional term directly changes the effective potential. For a pure SU{N) gauge 
theory, the effective potential Vg/f can be written as a character expansion of the form 

Veff = J2vnTrnP (2) 
where the sum is over all representations of zero A^-ality, i.e., invariant under Z(N). Terms of 



this form can be induced at one loop by certain topological excitations {2, 3, 4] as well as by 
particles in the adjoint representation. A one-loop calculation shows that the contribution 
to the effective potential of a heavy particle of mass M in the adjoint representation, either 
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boson or fermion can be approximated in 3 + 1 dimensions as 

'{2s + l)M^T^ 



-K2 (M/T) 



Tta (P) = -ThATvA (P) (3) 



where T is the temperature and 2s + 1 accounts for spin j^. The parameter Ha is positive in 
this case. The effect of such particles can be included at lowest order in Ha in the effective 
potential by the shift va —>■ va — THa- A positive value of Ha favors the Z(A^)-breaking 
deconfined phase. However, a term with Ha negative favors minimization of TtaP- Because 
TtaP = [FrpP] — 1, the minimization TtaP of implies TrpP = 0, a definining property of 
the confined phase. It is reasonable to expect that a sufficiently negative value of Ha might 
lead to a restoration of confinement at temperatures above the deconfinement temperature. 

We were motivated to look for this symmetry restoration by recent theoretical work on 
various aspects of the Polyakov loop effective potential. In certain supersymmetric gauge 
theories on B? ^ , Davies et al. 2|, yj have shown that finite temperature monopoles give 
rise to a Polyakov loop effective potential that has a Z(A^)-symmetric minimum for all values 
of the 5*^ circumference, and is therefore in a confined phase. These models do not precisely 
represent systems at finite temperature, because the supersymmetric partners of the gauge 
fields obey periodic boundary conditions. Comparable calculations in non-supersymmetric 
SU{N) gauge theories at finite temperature are much more difficult. In SU (2) gauge theory, 
Diakonov et al. [4] have calculated the contribution to Veff of finite-temperature instantons 
with non-trivial holonomy; such instantons have a color magnetic monopole content. Their 
work indicates an instability of the deconfined phase at sufficiently low temperature. In 
both of these examples, topological excitations give rise to a term in the effective potential 
corresponding to Ha negative. 

A positive value of Ha decreases the deconfinement temperature. For negative values of 
Ha, we have found new phases for both SU{3) and SU{A). In the case of SU{3), the new 
phase breaks Z{3) symmetry in an unfamiliar way, characterized by a negative value for the 
Polyakov loop in the fundamental representation (TrpP) < 0. In the case of SU{A), the 
global Z(4)symmetry is spontaneously broken to Z{2). The residual Z{2) symmetry ensures 
that for the fundamental representation (TrpP) = 0, but (TruP) 7^ for representations R 
that transform trivially under 2'(2), such as the 6 and the 10. 
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FIG. 1: SU{3) phase diagram in the /3 — plane. The dotted hne represents a plausible extrap- 
olation. 

SIMULATION RESULTS FOR SU{3) 



The lattice action we have studied for SU{3) and S'f/(4) is 



(4) 



where Sw is the Wilson action, defined conventionally as the sum over plaquettes. The 
sum in the second term is over all spatial sites, and naively Ha = hAd^- Most of our 
simulations were performed on 24^ x 4 lattices as reported here, but similar results were 
obtained for Nt = 2 and 6. The programs used for these simulation were developed us- 
ing the programming framework FermiQCD [6]. Because the augmented lattice action S 
depends quadratically on the time-like link variable Uq via the adjoint representation, the 
efficient heatbath methods developed for the standard lattice action cannot be used. We 
have used instead a recently developed SU{N) overrelaxation algorithm [7] combined with 
the Metropolis algorithm. The overrelaxation algorithm, which operates on the full SU (N) 
group rather than subgroups, proved to be fast and effective. Other algorithms which have 
been developed for fundamental plus adjoint actions could also be used s], Q]. A typical 
simulation on a 24^ x 4 lattice consisted of 10,000 equilibration sweeps followed by 60,000 
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FIG. 2: ^[/(S) Polyakov loop histogram at (3 = 6.5, Ha = -0.05. 
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FIG. 3: 5C/(3) Polyakov loop histogram at /3 = 6.5, //yi 
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FIG. 4: 5C/(3) Polyakov loop histogram at /? = 6.5, 



-0.08. 



sweeps during which 2000 measurements were performed. 

The approximate phase diagram for SU{3) is shown in Fig. [1] for A^^ = 4. The order 
parameter is TrpP, projected for each lattice field configuration onto the nearest Z{3) axis, 
with care to preserve the sign. There are three distinct phases: a deconfined phase where 
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FIG. 5: SU{3) Polyakov loop histogram at /? = 6.5, Ha = -0.1. 
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FIG. 6: 5*1^(3) Polyakov loop histogram at /? = 6.5, ff^ = -0.11. 
the projected expectation value satisfies {TrpP) > 0, a confined phase where (TrpP) = 0, 
and an intermediate phase with (TrpP) < 0, which we refer to as the skewed phase. The 
locations of the phase transitions were determined from the peaks of the adjoint Polyakov 
loop susceptibility, and checked against the histograms of the fundamental Polyakov loop. 
The dashed line in the phase diagram is an extrapolation; the phase transition between the 
skewed and confined phases is very difficult to resolve in this region. We will use the notation 
Hci for the values of Ha on the boundary between the deconfined and skewed phases, and 
Hc2 for the boundary between the skewed and confined phases. 

Figures [2] though [6] show histograms of the order parameter (TrpP) for various values 
of Ha a.t (3 = 6.5. At Ha = —0.06, there is clear evidence for the new intermediate phase 
where (TrpP) < 0. The skewed phase breaks Z{3) symmetry, as shown clearly by the 
histogram at Ha = —0.08, where all three possible skewed phases appear. The appearance 
of significant tunneling between the three phases on a 24^ x 4 lattice is an indication that 
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the transition from the skewed phase to the confined phase is very weak. The skewed phase 
differs from the deconfined phase not only in the orientation of the histograms, but also 
in the smaller magnitude of (TrpP) for the skewed phase. Near i^ci, the orientation of 
fiuctuations in histograms of the skewed phase is predominantly tangential, but becomes 
more radial as Hc2 is approached. The transition between the deconfined and skewed phase 
is clearly first-order, because the order parameter shows a marked jump when changing sign. 
The transition between the skewed phase and the confined phase is likely to be first order, 
because it is associated with the universality class of the three-dimensional Potts model 
and its generalizations via Svetitsky-Yaffe universality. However, in simulations (TrpP) 
shows a very small change at the skewed-confined transition, particularly near the apparent 
tricritical point. Empirically, for a given value of A^^^, the skewed phase shows up clearly 
only for Ns/Nt > 6. On a 12^ x 6 lattice, for example, the skewed phase always appears to 
coexist with either the deconfined phase or the confined phase. A detailed finite-size scaling 
analysis on very large lattices would be required to resolve the order of this transition with 
confidence. 

Figure [7] shows the projected value of (TrpP) for various values of Ha at /? = 6.5. The 
presence of three distinct phases is clear. The adjoint susceptibility xai for P = 6.5 is shown 
in Fig. [HI There is a clear peak between the deconfined and skewed phases, and a much 
smaller peak separating the skewed and confined phases. 



THEORY FOR 5C/(3) 

A simple theoretical approach based on the effective potential Veff for Polyakov loop 
eigenvalues reproduces the phase structure observed in simulations for SU{3) and SU{A). 
The effective potential has two parts. The first part is the one-loop expression for the free 
energy of gluons moving in a non-trivial, constant Polyakov loop background. The one-loop 



ree energy density was first evaluated by Gross, Pisarski, and Yaffe 



10|, and by N. Weiss 



llj . It is convenient to work in a gauge where Aq is a constant element of the SU{N) 
Lie algebra so that the background Polyakov loop is given simply by P = exp(i/3ylo)- The 
second contribution to the effective potential in our model is simply the term —hATTvAP 
that we have added to the gauge Lagrangian. At temperature T, our expression for Veff is 
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where the sum is over Matsubara frequencies c<j„ = 2imT. A useful form is 



N 



j,k=i 



T 



N 



- 1 



(6) 



where the angles 9j are the eigenvalues of PAq and|A6'jfc| is \9j — 9k\mod2TT. Thus 
is the sum of a one-loop term plus another term treated classically. 

The phase diagram is found by minimizing Vg/ / as a function of the Polyakov loop eigen- 
values. The two terms that make up Veff have identical local extrema, and the problem of 
minimizing Veff can be reduced to finding the minimum over this set. In the case of SU{3), 
it is sufficient to consider V^jj as TrpP varies along the real axis. In this case, the eigenval- 
ues of P may be taken to be the set {1, exp(20), exp(— ^0)}, and TrpP may be written as 
1 + 2cos(0). The effective potential is given by 



eff\ 



>T, h, 



Gtt' 



:(80^(0 - ny + <P\<P - 2ixY) - /iaT((1 + 2 cos(0))^ 



1 . 



(7) 



The extrema of V^ff occur at = 0, = 27r/3, and (p = ir. The values of TrpP for these 
values of are 3, 0, and —1, and we identify them with the deconfined, confined, and skewed 
phases, respectively. The set of eigenvalues {1, exp(27rz/3), exp(— 27rz/3)| is the unique set 



invariant under global Z{3) transformations 12l. Il3|. 



It is clear that the phase structure depends only on the dimensionless variable Ha/T^- 
As hj^ is lowered from zero, there is a first-order transition from the deconfined phase 
to the skewed phase. Setting the effective potential at = and = vr equal, we 
find that the transition from the deconfined phase to the skewed phase takes place at 
hci/T^ = — 7r^/48 ^ —0.206. As Ha decreases, another first-order transition, this time be- 
tween the skewed and confined phases, occurs at hc2/T^ = — 57r^/162 ^ —0.305 .We plot the 
potential as a function of Tr pP for values in the three phase in Figures [QlfTTl corresponding 
to Ha/T^ = 0, -0.24, -0.35. 

We cannot directly relate Ha and the corresponding lattice parameter Ha, because there 
is an unknown multiplicative renormalization relating the two. However, the ratio /ic2/^ci is 
approximately 1.48. If we assume that the relation of /i to if is approximately independent 
of h, we can compare with the results obtained from simulation. As shown in Fig. [H the 
ratios Hc2/Hci obtained vary from 1.27 at /5 = 6.2 to 1.44 aX (3 = 6.8, with a maximum 
value of 1.73 in between. 
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FIG. 9: Effective potential versus TrpP for deconfined phase aX Ha/T^ = 0- 
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FIG. 10: Effective potential versus TrpP for skewed pfiase at Ha/T^ 
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FIG. 11: Effective potential versus TrpP for confined phase at Ha/T^ = —0.35. 
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FIG. 12: Theoretical prediction for pressure normalized to black body pressure pressure as a 
function of Ha- 

As noted previously, our simulations show a pronounced asymmetry in the skewed phase 
between the fluctuations of the imaginary and the real parts of TrpP- Fluctuations in the 
projected imaginary part are associated with motion in the Ag direction, while fluctuations 
in the projected real part are due to motion in both the As and A3 directions. It is thus 
interesting that in the skewed phase, theory predicts an asymmetry in the screening masses 
obtained from small fluctuations in the eigenvalues of P. This is quite different from the 
behavior in the conflned and deconflned phases, where theory predicts no asymmetry. We 
have 

ma^ / 1 + 2hA/T^ 
ms 



(8) 



-1 - 6hA/T^ 

This ratio varies from 1.59 at hd to 0.69 at hc2- This is on the order of the variation seen 
in the fluctuations of the real and imaginary parts of TrpP, and probably accounts for the 
behavior seen in the histograms. This prediction for the mass ratio can be checked more 
directly by comparing the masses obtained from the correlation functions of the real and 
imaginary parts of the projected Polyakov loop in the skewed phase. 

The pressure can be calculated from simulations along a path of constant (3, using 

r-2 



Y^-Y^-N'n dHAiTTAP) 



(9) 
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A detailed comparison of the pressure for all values of Ha would require knowledge of the 
relation between Ha and Ha- However, it is relatively simple to compare the change in the 
pressure from Ha = Ha = to the deconfined-skewed phase boundary as well as the change 
in pressure across the skewed phase. Using Veff, we find that the predicted change in p/T^ 
from hci to is 7r^/6 ~ 1.64; from hc2 to the net change is 0. For comparison, the 
corresponding results from simulations at /? = 6.5 are 1.64 ± 0.03 and P = —0.18 ± 0.07. In 
each case, the error is completely dominated by systematic error due to uncertainty in the 
location of the critical values of Ha, with statisitical error at least an order of magnitude 
smaller. 

SIMULATION RESULTS FOR SU{A) 

We have also simulated SU{A) lattice gauge theories, again primarily on 24^ x 4 lattices. 
As in the case of SU{3), we find a new phase in the region Ha < 0, but the nature of 
the new phase is completely different. In this new, partially confined phase, global 2'(4) 
symmetry is spontaneously broken to Z(2). In this phase, particles in the fundamental 
representation ("5'?7(4) quarks") are still confined, but bound states of two such particles 
("S'f/(4) diquarks") are not. Each irreducible representation of SU{N) has an A^-ality: if 
z G Z{N), P —>■ zP induces a change Trj^P — > z^Tr^P, where k is the A^-ality of the 
representation R. The characteristic feature of the partially confined phase in SU{4) is that 
the expected value of Polyakov loops in = 1 representations is zero, but not in A; = 2 
representations such clS clS the 6 and the 10. 

The breaking of Z{4) down to Z{2) for sufficiently negative Ha is manifest in histograms 
of the Polyakov loop in the fundamental representation as a clustering of data around either 
the X OT y axis, but not both, as shown in Figures fTSlfTBl The Z{2) character of this new phase 
is very clearly shown in Figure [TTl which shows the behavior of the real and imaginary part 
of the Polyakov loop versus Monte Carlo time for one long run with 20,000 measurements. 
As the figure reveals, there are significant fluctuations in either the real or the imaginary 
part, but not both simultaneously, characterisitic of Z{A) breaking to Z{2). In this phase, 
the expectation value of TrpP^ is non-zero, being positive when the fluctuations in TrpP 
are along the real axis, and negative when TrpP fluctuates along the imaginary axis. 

As Ha becomes more negative, the histograms show decreasing amplitude in the fluctu- 
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FIG. 13: S'C/(4) Polyakov loop histogram at /? = 11.1, Ha = -0.1. 
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FIG. 14: SC/(4) Polyakov loop histogram at /3 = 11.1, Ha = -0.11. 



ations of TrpP- It is possible that there is a second phase transition from the Z{2) phase to 
the confined phase as Ha becomes more negative, but we have not found direct evidence for 
this. As we discuss below, our simple theoretical model does not predict a second transition 
for this theory, at least not at high temperatures where it is valid. 




FIG. 15: SU{A) Polyakov loop histogram at /3 = 11.1, Ha = -0.12. 
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FIG. 16: 5[/(4) Polyakov loop histogram at /? = 11.1, Ha = -0.125. 
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FIG. 17: Real and imaginary parts of S'f7(4) Polyakov loop versus Monte Carlo time at /? = 11.1, 
Ha = -0.11. 

THEORY FOR 5f/(4) 



We have examined within our simple theoretical model the possible occurrence of four 
different phases in SU{4): the confined phase, which has full Z{4) symmetry; the deconfined 
phase; a partially-confined, Z(2)-invariant phase; and a skewed phase similar to the skewed 
phase of SU{3). Only the deconfined phase and the Z{2) phase are predicted by our simple 
theoretical model. 

The properties of the Z(2)-invariant phase may be understood by considering the one- 
parameter class of eigenvalues invariant under Z{2); the eigenvalues in this class may be 
written as {9,7t — 6,tt + 6, 2tt — 6}, and the corresponding Polyakov loops have the form 
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diag [e*^, — e — e'^, e . The one-loop effective potential as a function of 9 becomes 

Veff = + ^ [(^ - 2^)' + 2^)' + + (27r - 2^)' (2^)'] + HaT (10) 

which has its minimum within this class at 6* = 0. The confined phase, which has Z{A) 
symmetry, is realized at ^ = 7r/4 but is never the minimum of K//. This behavior is easy 
to understand: both the confined and Z(2)-invariant phases have the same dependence on 
Hai so the stable phase is the one that minimizes the contribution of the gauge bosons. The 
deconfined phase does not fall into the Z(2)-invariant class: with all eigenvalues set to 0, 
the value of the effective potential in the deconfined phase is 

Vd = — - IhhAT (11) 

There is a first-order transtion between the deconfined and Z(2)-invariant phases at Ha/T^ — 
-7rV48 ~ -0.205617. The value of A {p/T^) between Ha = ^ and the critical point is 
7r^/3 ~ 3.289. The value we obtained from simulations at /3 = 11.0 was 2.21 ± 0.07, where 
again the systematic error dominates, due to uncertainty in the location of the transition. 

In order to realize the confined phase, it may be necessary to add an additional term 
proportional to TtaP"^ = TrpP'^TrFP^^ - 1 in order to force both TrpP and Tr^P^ to 
zero, but this has not yet been checked in simulations. Of course, there must be a line of 
transitions in the (3 — Ha plane separating the Z{2) phase from the low-temperature confined 
phase with Z{4) symmetry. The transition could be either first or second order. We have 
not yet mapped out this phase boundary via simulation. As previously noted, our simple 
theoretical model does not include a mechanism for this transition. 

CONCLUSIONS 

We have considerable evidence, from lattice simulation and from theory, for the existence 
of new phases of finite temperature gauge theories, and for the restoration of the confined 
phase at high temperatures when extra, Z(A^)-invariant, Polyakov loop terms are added to 
the gauge action. In SU{3), a novel skewed phase was found, and in S'[/(4), we found a 
phase where Z(4) is spontaneously broken to Z{2). In the general case of SU{N), there is 
good reason to expect a very rich phase structure may exist. 
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A simple theoretical model based on perturbation theory at high temperatures has proven 
surprisingly accurate in predicting the observed phase structure and thermodynamics. Al- 
though successful, the model has significant shortcomings. Fluctuations in Aq are not con- 
sidered, nor is the renormalization of h^- Most importantly, our simple model does not 
include in Vgff whatever mechanism is responsible for confinement at low temperatures 
when Ha = 0. It is therefore invalid at low temperatures. Nevertheless, theory and simula- 
tion are in reasonable agreement on a wide range of properties. Our model can also make 
predictions for string tensions and 't Hooft loop surface tensions, and these predictions can 
be checked in lattice simulations. The introduction of Ha as an extra parameter also affects 
the action of calorons, topologically stable solutions of the classical field equations, and thus 
may offer rich possibilities for explorations in instanton physics. 

The interpretation of these additional phases of finite temperature gauge theories is, to 
a degree, associated with the issue of a physical implementation of a negative value for 
Ha- The existence of a partially confining, Z(2)-invariant phase in SU{4) might have been 
expected 121], and the interpretation of the order parameters is clear. The interpretation of 
the skewed phase in SU (3) is less certain. As in the deconfined phase, the global symmetry 
of the Polyakov loop is lost in the skewed phase. Our theoretical analysis indicates that, 
on average, two of the three Polyakov loop eigenvalues are degenerate, suggesting a possible 
interpretation of the skewed phase as some form of SU{2) x f/(l) Higgs phase. 

The issues underlying the interpretation of parameters and phases is connected with the 
association of finite temperature gauge theories with universality classes of spin systems ^ . 
It has always been assumed implicitly that the mapping from gauge theories to spin systems 
is into but perhaps not onto. There are phases of SU{N) and Z{N) spin systems which 
are not easily obtainable from physical finite temperature gauge theories. For example, the 
antiferromagnetic phase of a spin system can be obtained from the strong-coupling effective 
action of a lattice gauge theory with < and Nt odd, a construction with no obvious 
continuum limit. However, phases can often be reached in different ways in the space of 
parameters. The skewed phase we have found in SU{3) gauge theory is very similar to 
the anti-center phase found in SU{3) spin systems by Wozar et al. ij]. Although the 
term in the spin Hamiltonian that produces the anti-center phase is associated with the 15 
representation rather than the adjoint term we have used, we are confident that the two 
phases will prove to be related. At this time, it is simply unclear what physical principles. 
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if any, limit the map between spin systems and gauge theories. 

We beheve that the abihty to create new phases in a controlled way may become an 
important tool in understanding the properties of finite temperature gauge theories. For 
example, simulations indicate that the confined phase obtained in SU{3) at high T with 
< is connected to the conventional confined phase at low T with = 0. The 
possibility of a confined phase in a region where perturbation theory is valid is by itself 
enormously interesting. As larger gauge groups are considered, the number of possible new 
phases increases. For example, in SU{6), we can consider partial breaking of Z{6) to either 
Z{2) or Z{3). At high temperatures, analytic calculations of both string tensions and 't Hooft 
loop surface tensions can be carried out in these different phases for potential comparison 



with simulation 
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